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ABSTRACT 

(Distribution  Limitation  Statement  No.  2) 

Using  Fresnel  diffraction  integrals,  calculations  have  been  made  of  the 
irradiance  and  power  distributions  in  the  vicinity  of  the  focus  for  Gaussian 
beams  focused  through  annular  apertures.  The  effect  of  shifting  the  peak  of 
the  Gaussian  aperture  distribution  from  the  aperture  center  was  investigated 
and  was  found  to  be  a  small  effect  in  some  special  cases.  Gaussian  and 
sinusoidal  phase  distortions  of  the  aperture  field  were  included  in  the  calcu¬ 
lations.  It  was  found  that  these  kinds  of  phase  distortions  generally  reduce 
the  focal  plane  irradiances  in  sometimes  unpredictable  ways.  The  numerical 
examples  chosen  for  illustration  should  be  useful  in  interpreting  future 
experiments  to  be  conducted  by  the  Air  Force  Weapons  Laboratory  at  the  Sandia 
Optical  Range. 
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SECTION  I 
INTRODUCTION 

In  this  report,  we  consider  some  of  the  features  of  focusing  10.6  micron 
laser  beams  at  kilometer  ranges.  To  use  a  simple  analytical  approach  we  assume 
that  the  focusing  device  is  an  ideal  on-axis  Cassegrain  telescope  using  spherical 
mirrors  and  that  the  field  amplitude  is  Gaussian  in  magnitude  over  the  final 
transmitting  mirror  of  the  telescope.  Various  simple  kinds  of  phase  distortion 
are  introduced  and  the  subsequent  changes  in  the  nature  of  the  focused  beam  are 
discussed.  Numerical  illustrations  have  been  prepared  using  specific  parameter 
values  rather  than  generalized  dimensionless  coordinates.  The  beam  propagation 
is  assumed  to  be  diffraction  limited;  all  other  possible  propagation  factors, 
such  as  medium  turbulence,  are  Ignored  in  this  work. 
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SECTION  II 

DIFFRACTION  FORMULATION 

To  set  up  the  problem,  we  define  an  (r,  6,  z)  cylindrical  coordinate  system. 
The  final  transmitting  mirror  of  the  telescope  is  assumed  to  lie  in  the  z  =  0 
plane  while  the  output  beam  travels  in  the  positize  z  direction.  We  assume  that 
the  focusing  problem  can  be  handled  by  postulating  a  convergent  Gaussian  beam 
diffracted  by  an  annular  aperture  in  the  z  =  0  plane.  The  transmitting  portion 
of  the  annular  aperture  occupies  the  region  a  £  r  £  b.  Physically,  this  can  be 
interpreted  as  a  collimated  beam  of  radius  a,  entering  the  telescope  and  emerging 
as  an  annular  focused  beam  of  inside  radius  a  and  outside  radius  b,  at  the  final 
mirror  of  the  telescope.  The  diffraction  calculations  are  applicable  to  the 
region  z  >  0  and  are  subject  to  the  usual  restrictions  governing  the  use  of 
Fresnel  diffraction  integrals  in  the  scalar  wave  approximation.  The  case  of  the 
off-axis  Cassegrain  telescope  can  be  considered  by  simply  setting  a  =  0  in  which 
case  a  is  not  interpreted  as  the  radius  of  the  input  beam  to  the  telescope. 

The  aperture  field  is  taken  as 

v(r)  =  vQ  exp[-(r/w)2-iTTr2/\R-i<J>]  (1) 

In  the  sense  of  geometrical  optics,  the  telescope  is  focused  at  a  point  on  the 
z-axis  a  distance  R  from  the  z  =  0  plane,  hence,  the  plane  z  =  R  is  called  the 
focal  plane.  The  focal  range  R  is  taken  to  be  positive  which  implies  a  time 
variation  given  by  exp[-iwt],  In  (1),  w  is  the  beam  spot  size  at  the  telescope 
transmitter  mirror,  X  is  the  vacuum  wavelength,  and  4  is  a  phase  distortion  which 
we  shall  restrict  to  be  a  function  of  r  only.  For  most  of  the  calculations 
described  herein,  the  peak  of  the  Gaussian  distribution  is  on  the  z-axis  as 
indicated  by  equation  (1).  In  one  later  section,  however.  Equation  (1)  will  be 
modified  so  that  the  peak  of  the  Gaussian  is  shifted  to  the  point  (r,  0)  = 

(Ax,  o). 

For  convenience  in  later  calculations,  we  shall  assume  that  v  is  a  real 

’  o 

constant  and  is  normallized  so  that  total  power  PQ  is  always  transmitted  by  the 
telescope.  In  the  scalar  wave  approximation,  we  can  then  always  determine  vQ 
by  the  relation. 
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D 

PG  =  2tt  J  dr  r  |v(r)  | 


(2) 


In  all  of  the  figures  to  follow,  we  take  PQ  =  150  watts.  When  v(r)  is  given  bv 
equation  (1)  we  have 


vo  =  (2P0/ttw2  j/ [exp(-2  (a/w)2)-exp(-2(b/w2)] 


(3) 


Using  the  Fresnel  diffraction  integral,  the  scalar  wave  complex  amplitude  of 
the  diffracted  beam  is 

b 

u(r,  z)  =  (2ir/Az)  J  dp  p  v(p)  J0(2irpr/Az)exp[iirp2/Az]  (4) 


where  we  have  omitted  a  phase  factor  on  the  right  hand  side  of 


-i  exp[i27tz/A+iTTr2/Az] 


The  irradiance  of  the  diffracted  beam  is  taken  simply  as  I(r,  z)  =  |u(r,  z) |  . 
In  any  plane  z  =  constant,  the  power  passing  through  a  circle  of  radius  r  centered 
on  the  z-axis  given  by 

_r 

.  2 

(5) 


r 

P (r ,  z)  =  2tt  J  dp  p  | u(p ,  z)  | 
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SECTION  III 
ON-AXIS  IRRADIANCE 


We  now  assume  <t>  =  0 
restrictions  1(0, z)  can 


in  equation  (1)  and  r  =  0  in  equation  (4). 
be  written  as 


With  these 


where 


1(0, z)  =  v 


?  NUM 
o  DEN 


NUM  =  [exp (- (a/w) 2 )  -  exp(-b/w)2)] 


20  2 


+  4  exp[-(a2+b2)/w2]  sin2  [it  (b2-a2)  (z-R)/2AzrJ 


DEN  =  (Az/ttw2)2  +  (z-R)  2/R2 


(6) 


It  is  commonly  assumed  that  l(0,z)  is  symmetrical  about  the  point  z  =  R  but 
inspection  of  equation  (6)  reveals  that  this  is  not  generally  the  case. 

There  is  an  optimization  problem  that  one  can  readily  solve  using  equation 
(6).  We  first  consider  that  z  is  fixed  at  a  constant  value  zq  and  that  R  is 
variable.  In  this  case  I(0,zQ)  is  regarded  as  a  function  of  R  and  there  exists 
an  optimum  value  of  R  which  maximizes  I(0,zo).  This  optimum  value  is  simply 
R  =  zQ.  Physically,  this  means  that  an  on-axis  detector  located  at  an  arbitrary 
fixed  location  zQ  will  receive  the  largest  signal  when  the  telescope  is  focused 
such  that  R  =  zQ.  This  is  true  regardless  of  the  values  of  a,  b  and  w.  If  now 
the  telescope  focus  is  fixed  at  R  =  zQ  and  the  on-axis  detector  is  physically 
moved  along  the  z-axis  toward  the  telescope,  it  is  generally  found  that  the 
detector  signal  will  first  increase,  then  pass  through  a  maximum  and  then  decrease. 

Figures  1,  2,  and  3  illustrate  some  of  the  predictions  of  equation  (6).  All 
three  figures  are  calculated  for  a  wavelength  of  10.6  microns.  In  figures  2  and 
3,  the  totally  unobscured  or  complete  Gaussian  beam  is  shown  as  the  solid  line 
for  comparison.  The  focal  plane  is  located  by  the  arrow  on  the  abscissa.  It  is 
interesting  to  note  that  the  focal  plane  irradiance  for  each  case  represented  by 
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a  broken  line  is  approximately  one-half  of  the  focal  plane  irradiance  for  the 
free  Gaussian.  Both  on-axis  and  off-axis  Cassegrain  telescopes  were  used  to 
penetrate  these  curves.  The  figures  verify  that,  for  fixed  a,  b,  P  ,  and  A,  the 
focal  point  irradiance  is  maximized  for  w  i.e. ,  the  case  of  uniform  illumi¬ 

nation.  An  analytical  proof  of  this  maximization  is  as  follows. 

We  evaluate  the  on-axis  irradiance  at  the  focal  point,  define  1^  =  1(0, R), 
and  obtain 

2ttw2P  exp  [(b2-a2  ) /w23 -1 

t  =  - °  .  - - -  (7) 

(AR)2  exp  [(b  -a2)/w23+l 

Taking  the  partial  derivative  of  1^  with  respect  to  w,  we  can  obtain 


31,  WP  I i  [(b2-a2)/w2]n[2n-2n]/n! 

J-- - 2-  '  2=% - ry—  (8) 

(AR)2  ^exp  [(b2-a2) /w23+l^ 

Now,  since  2n  >  2n  for  n  >  3,  we  see  that  3I^/9w  is  always  positive  and  therefore 
1^  increases  as  w  increases.  The  maximum  value  of  1^.  occurs  in  the  limit  w  -►  <*> 
and  is  given  simply  by 


lim  If  =  ir(b2-a2)Po/(AR)2  (9) 


Note  that  we  have  restricted  the  above  analysis  to  the  case  where  the  total  beam 
power  transmitted  by  the  telescope  is  fixed  at  the  constant  value  P0. 
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Figure  1.  I(o,z)  as  a  Function  of  z.  Calculated  from  Equation  (6) 

(The  ordinate  has  units  of  watts/cm2  when  PQ  =  150  watts.  The  focal  plane  is 
located  by  the  arrow  on  the  abscissa.) 
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Figure  2.  1(0, z)  as  a  Function  of  z,  Calculated  from  Equation  (6) 

(The  ordinate  has  units  of  watts/cm2  when  the  total  beam  power  P0  =  150  watts. 
The  wavelength  is  10.6  microns.) 
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of  z,  Calculated  from  Equation  (6) 
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SECTION  IV 
TRANSVERSE  PROFILES 


We  shall  now  restrict  our  attention  to  the  focal  plane  and  examine  the 
transverse  irradiance  and  power  profiles  for  ideal  beams  with  4>  =  0.  Under 
these  restrictions  equation  (4)  reduces  to 

b 

u(r,R)  =  (2irvo/XR)  J  dp  p  JQ(2TTrp/XR)exp  [-(p/w)"J  (10) 

a 

There  does  not  appear  to  be  a  simple  way  to  evaluate  equation  (10)  in  terms  of 
elementary  tabulated  functions,  hence,  numberical  illustrations  were  prepared 
by  computer  integration  of  equation  (10).  Transverse  irradiance  profiles 
I(r,R)  =  |u(r,R)|2  are  shown  in  figures  4  and  5.  The  free  Gaussian  is  again 
shown  for  comparison  purposes.  The  curve  for  a*o  and  w  =  “  is  the  familiar  Airy 
disk  pattern  whose  first  dark  ring  (zero)  occurs  for  r  -  0.61  XR/b.  The  values 
chosen  for  a,  b,  w,  X,  R  and  Pq  are  the  same  as  for  figures  2  and  3.  The  evalua¬ 
tion  of  vQ  is  given  by  equation  (3). 

The  transverse  power  profiles  (normalized  to  PQ)  are  given  by 

P(r)/PQ  -  (2*/P0)  f  dp  p  | u (p , R) | 2  (11) 

o 

and  are  shown  in  figures  6  and  7. 
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10 


r  (cm) 

Figure  4.  Transverse  Irradiance  Profile  I(r,R)  as  a  Function  of  r 

(The  ordinate  has  units  of  watts/cm2  when  the  beam  power  P  =  150  watts. 
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Figure  5.  Transverse  Irradiance  Profile  I(r,R)  as  a  Function  of 
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Figure  6.  Normalized  Focal  Plane  Power  P(r)/PQ  (in  percent)  as  a  Function  of  Focal  Plane  Radius 

(The  effects  of  the  central  obscuration  in  the  on-axis  Cassegrain  telescope  are  most 
noticeable  in  the  60  to  90  percent  range.) 


AFWL-TR- 70-146 


r  (cm) 

Figure  7.  Normalized  Focal  Plane  Power  P(r)/P0  (in  percent)  as  a  Function  of  Focal  Plane  Radiu 
(The  effects  of  the  central  obscuration  in  the  on— axis  Cassegrain  telescope  are  most 
noticeable  in  the  75  to  90  percent  range.) 
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SECTION  V 

GAUSSIAN  PHASE  DISTORTION 

In  this  section  we  consider  the  case  when  the  phase  distortion  in  equation 
(I)  is  given  by 

4>  =  (2rA/X)exp [-2 (r/w)2}  (12) 

The  motivation  for  assuming  $  of  the  form  equation  (12)  is  to  provide  a  simple, 
though  perhaps  crude,  model  of  the  effects  owing  to  beam  heating  of  mirror 
surfaces.  Every  mirror  in  the  optical  train  prior  to  the  telescope  and  even 
including  the  telescope  mirrors  absorbs  a  small  fraction  of  the  beam  power. 

This  small  fraction  is  not  insignificant,  however,  because  the  absorbed  energy 
is  converted  to  heat  near  the  surface  of  the  mirror  and  this  heat  forces  dimen¬ 
sional  changes  in  the  mirror  surface.  The  simplest  possible  model  is  to  postu¬ 
late  that  each  mirror  expands  in  a  direction  normal  to  its  reflecting  surface 
by  an  amount  proportional  to  the  local  beam  irradiance.  Thus  a  normally  plane 
mirror  surface  will  become  a  convex  Gaussian  surface  when  subjected  to  a  Gaussian 
local  irradiance  distribution. 

By  neglecting  near  field  diffraction  between  mirrors,  one  then  arrives  at  the 
Gaussian  phase  term  as  given  in  equation  (12).  Such  a  phase  distortion  can  be 
visualized  as  consisting  of  two  significant  components:  a  quadratic  term  in  r 
plus  a  more  complex  function.  Thus  we  could  write  equation  (12)  as 

*  *  (2ttA/X)-|  Trr2 /XR'}+ |  (2TtA/X)exp  [-2  (r/w)2]- (2TTA/X)+(rr2/XR')}  (13) 

where  R'  is  not  specified  as  yet.  In  equation  (13),  the  first  term  in  braces  |J 
can  be  interpreted  as  a  distortion  component  that  only  changes  the  focal  point 
(or,  more  properly,  image  point)  of  the  telescopic  system.  If  w  is  comparable 
to  b,  it  is  found  numerically  that  there  exists  an  R'  such  that  the  second  term 
in  braces  in  equation  (13)  has  only  a  small  effect  on  the  diffraction 
calculations . 
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Since  the  Gaussian  distortion  model  is  only  approximate  and  perhaps  somewhat 
naive,  no  extensive  effort  has  been  expended  to  determine  validated  analytical 
solutions  for  R'.  Possibly  a  simple  least  squares  fit  of  4-2ttA/A  to  the  para¬ 
boloid  -Trr2/AR^  would  suffice.  In  this  case,  R'  would  then  be  given  by 


R^  =  (it/ A )  J'  dr  r5  /  J'  dr  r 3  (2ttA/A— 4?) 


(14) 


The  upper  limit  of  the  integrals  in  (14)  would  be  b  while  some  decision  would 
have  to  be  made  whether  to  make  the  lower  limit  zero  or  a.  The  new  focal  point 
(or  image  point)  would  then  be  given  by  R"  where 

1/R"  =  1/R  -  1/R"  (15) 

and  R"  is  considered  positive  for  positive  A. 

An  alternative,  and  simpler,  technique  for  determining  R"  is  simply  to  expand 
4>  in  an  infinite  series,  yielding 


4>  =  (2-rrA/A)  [l-2(r/w)2  +  higher  order  terms]  •  (16) 

The  quadratic  term  identifies  R"  as 

R'  ■=  w2/4A  (17) 

Equation  (17)  has  the  advantage  of  being  simpler  than  equation  (14)  and  for 
w  i  b  is  probably  adequate.  The  effective  focal  point  is  given  simply  by 

1/R"  -  l/R-4A/w2  (18) 

Calculations  have  been  made  of  the  on-axis  irradiance  assuming  various 
numerical  values  for  A/A  in  equation  (12).  The  diffraction  integral  was  evaluated 
by  using  a  series  expansion  for  v(r)  rather  than  direct  computer  integration.  The 
diffracted  field,  evaluated  on  the  z-axis,  is  given  by 
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u(o,z)  =  (2nvo/Az)£]  [(— i 2tt A / A )n  /  n!] 

n=o 


r 

”1 

f  dp  p  exp  p 

2  |  -  (2n+l)w2+iiT  (R-z)  /ARz  | 

(19) 


After  evaluating  the  integral,  u(o,z)  becomes 


u(o,z)  =  (ttvq/Az)  X]  [(-i2TTA/A)n/n!  cnJ  [exp (-a2cn)-exp (-b2cn)]  (20) 

n=o 

where 

c  =  (2n+l)/w2  -Tii  (R-z)/ARz  (21) 

n 

Equations  (3),  (20),  and  (21)  were  used  to  prepare  figures  8,  9  and  10. 

2 

In  figure  8,  we  show  the  on-axis  irradiance  I(o,z)  =  |u(o,z) |  as  a  function 
of  distance  from  the  telescope  with  the  telescope  focus  fixed  at  1.7  km.  As  A/A 
increases  the  irradiance  at  z  =  1.7  km  decreases.  If  a  10.6  micron  beam  is 
turned  on  at  time  t  =  o  and  A  depends  on  time  thereafter  in  a  known  way,  these 
curves  could  possibly  be  useful  for  depicting  the  time  evolution  of  the  on-axis 
irradiance. 

Figure  9  shows  that  a  large  fraction  of  the  distortion  (Eq.  (12))  contributes 
to  an  effective  change  in  focus  of  the  telescopic  system  and  indicates  that 
equation  (17)  is  a  reasonable  approximation  for  w  >  b.  The  same  conclusion  is 
drawn  from  figure  10,  which  is  an  alternative  description  of  the  predictions  of 
equation  (20) . 

For  each  value  of  A/A  in  figure  10,  the  telescope  focus  R  was  adjusted  to 
locate  the  peak  of  the  on-axis  distribution  at  the  axial  point  z  =  1.7  km,  which, 
it  should  be  noted,  is  not  the  same  as  maximizing  the  irradiance  at  the  location 
z  =  1.7  km. 
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Figure  8.  On-Axis  Irradiance  1(0, z)  ■  |u(0,z)|  as  a 
Function  of  z,  Calculated  from  Equation  (20) 

(The  integer  n  labels  each  curve  and  also  specifies  the  value  of  A/A  associated 
with  that  curve.) 
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Figure  9.  On-Axis  Irradiance  at  a  Fixed  Distance  z  -  1.7  km  as  a 
Function  of  Telescope  Setting  R  for  Various  Values  of  A/A 


(The  integer  n  =  1,  2,  9  which  identifies  each  curve  specifies  the  value  of 

A/A  =  (n-l)/8. ) 
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SECTION  VI 

SINUSOIDAL  PHASE  DISTORTION 


In  this  section,  we  consider  the  case  when  ♦  is  given  by 

$  “  (2irm/A)cos  (2rr/T)  (22) 

where  m  specifies  the  peak  value  and  T  specifies  the  period  of  the  sinusoidal 
phase  distortion.  One  of  the  reasons  for  choosing  to  examine  equation  (22)  is 
that  the  quadratic  (change  of  focus)  component  should  be  small  for  T  <  b,  which 
allows  one  to  then  directly  investigate  the  degrading  effects  of  the  sinusoidal 
phase  distortion  on  the  focal  plane  irradiance  and  power  distributions.  In  an 
approximate  way,  equation  (22)  does  have  connection  with  physical  reality.  For 
example,  there  exist  mirrors  that  have  surface  deformation  that  depend  on  radius 
in  a  nearly  periodic  manner  while  possessing  azimuthal  symmetry. 

When  T  >>  b,  we  can  see  immediately  that  equation  (22)  represents  mainly  a 
change  of  focus  with  R'  found  from  an  expansion  of  equation  (22)  : 


-trr2/AR'  *  -  (2ttiii/A)  (1/2)  (2irr/T)2 


to  yield 


R'  -  T2/47i2m 


(23) 


In  this  section,  however,  we  are  mainly  interested  in  the  case  T  <  b.  For 
this  case  we  shall  assume  that  the  beam  induced  change  of  focus  is  negligible  so 
that  we  can  go  directly  to  the  telescope  focal  plane  R  ■  z  to  investigate  the 
transverse  distributions. 

To  portray  some  of  the  effects  of  the  cosinusoidal  phase  distortion 
(Eq.  (22))  we  proceed  directly  with  numerical  illustrations.  In  figure  11  we 
show  the  irradiance  at  the  focal  point  as  a  function  of  peak  distortion.  It  was 
found  that  the  focal  point  irradiance  was  independent  of  T  for  the  range  1  cm  i 
T  i  15  cm.  Note  that  the  focal  point  irradiance  is  nearly  zero  for  m/A  -  3/8  and 
T  -  1  cm.  The  entire  transverse  irradiance  and  power  distributions  are  reduced 
in  correspondence  to  the  focal  point  irradiance  for  this  case. 
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2 

Figure  11.  Focal  Point  Irradiance  1(0, R)  -  |u(0,R)|  as  a  Function  of 
Peak  Distortion  m/A  with  Period  T  as  the  Parameter  that 
Changes  from  Curve  to  Curve 

Figure  12  shows  the  transverse  irradiance  and  power  distributions  in  the 
focal  plane  for  a  period  T  =  1  cm.  Note  that  the  numbers  for  the  m/A  =  3/8  case 
are  too  small  to  show  up  on  the  scale  of  these  graphs.  The  power  distribution 
for  the  m/A  =  3/8  case  is  partially  shown  on  an  extended  scale  in  figure  13. 

In  figures  14  and  15,  we  show  irradiance  and  power  distributions  for  the  some¬ 
what  larger  spatial  periods  of  T  ■  15  cm  and  T  -  20  cm. 
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Figure  12.  Focal  Plane  Irradiance  and  Power  Profiles 
(The  period  T  =  1  cm  for  all  curves.) 
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Figure  14.  Focal  Plane  Irradiance  and  Power  Profiles 
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SECTION  VII 

DISPLACED  GAUSSIAN  PROFILES 

In  this  section,  we  assume  an  ideal  (<J>  =  0)  but  displaced  Gaussian  aperture 
distribution.  The  aperture  field  is  displaced  by  an  amount  Ax  along  the  x-axis 
and,  in  cylindrical  coordinates  (r,0),  is  given  by 

v(r,0)  =  vQexp  |  (rcos0-Ax)2+(rsin0)2  [ /w2-iirr2/AR]  (24) 

We  still  normalize  vQ  to  a  beam  power  PQ  so  that  vQ  is  found  from 


o  a 


(25) 


The  diffracted  wave  amplitude  is  defined  as  u(r,0,z).  We  shall  confine  our 
attention  to  the  focal  plane  z  *  R  so  that  u(r,0,R)  can  be  calculated  by  a 
Fraunhofer  integral,  as  follows 


2tt  b 

u(r,0,R)  -  (AR )~l  J  d<j>  J'  dp  p  |v(p,<J>)|  exp[-2irirp  cos(0-<f>)/AR]  (26) 


where  we  have  again  ignored  purely  phase  factors  on  the  right  hand  side  of 
equation  (26). 

Figure  16  shows  the  results  of  some  computer  integrations  of  equation  (25) 
and  equation  (26).  The  irradiance  profile  is  taken  along  the  x-axis  in  the 
focal  plane  where  the  displacement  Ax  has  the  greatest  effect. 


26 


=  25  cm 
=  1.7  km 
—  AX  * Ocm 


AFWL-TR- 70-146 


w 

x 

< 

X 

.23 

4J 

4J 

00  C 

c  <u 
o  a 

dj 

td  u 

m 

SI  ?— i 

—  n. 
<— ■>  to 
ci  -h 
•>  Q 
O 
*  T) 
M  i— I 
(u 
3  -H 
- U* 

II  0) 
M 

r-s  3 
PC  w 
«  M 
O  <U 
»  O. 

M  C 


O 

to 

tl)  to 

r-l  (U 

•rl  3 
IM  r-l 

O  m 

n  > 

C4 


Q) 

o 

e 

td 

•rl 

•o  > 

td 


to 

3 

O 

•H 

U 

td 


n 

u 


a) 

c 

td 
i — i 

pl. 


td 

u 

o 


SO 


<U 

Pj 

3 

00 

•H 


27 


AFWL-TR-70-146 


SECTION  VIII 
CONCLUSIONS  AND  SUMMARY 


This  report  contains  some  specific  calculations  of  Fresnel  diffraction 
integrals  using  parameter  values  that  should  be  useful  in  predicting  the 
diffraction-limited  propagation  of  Gaussian  laser  beams  launched  over  long 
ranges.  Both  ideal  and  phase  aberrated  aperture  distributions  have  been  con¬ 
sidered.  It  was  found  that  Gaussian  phase  distortions  of  small  value  generally 
result  in  an  effective  change  of  focus  in  the  telescopic  system  while  cosinu¬ 
soidal  phase  distortions  of  small  period  change  the  focal  plane  distributions 
in  generally  unpredictable  ways. 

The  on-axis  irradiance  profiles  show  that,  for  long  focal  length  infrared 
systems,  these  profiles  are  not  symmetrical  about  the  focal  point.  The  possi¬ 
bilities  for  asymmetries  of  this  nature  have  apparently  been  overlooked  in 
advanced  texts  on  optics.  For  example,  in  M.  Born  and  E.  Wolf,  Principles  of 
Optics ,  (1st  edition,  1959),  page  439,  we  find  the  statement  that,  "in  the 
neighborhood  of  the  focus  the  intensity  distribution  is  symmetrical  about  the 
geometrical  focal  plane." 
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